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ABSTRACT 

The photoevaporation of circumstellar disks is a powerful process in the disk dissipation at the origin of the Orion 
proplyds. This Letter reports the first detection of a photoevaporating disk in the final but long-lasting phase of its 
evolution. The disk is associated to a low-mass T Tauri member of the a Ori Cluster. It is characterized by a very 
low (if any) accretion rate and by a tenuous (Mioss^ 10~^ Mo/yr) photoevaporation wind, which is unambiguously 
detected in the optical spectrum of the object. The wind emits strong forbidden lines of [SII] and [Nil] because the 
low-mass star is close to a powerful source of ionizing photons, the 09.5 star a Ori. 
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^ 1. Introduction 



The evolution of circumstellar disks surrounding low-mass 
T Tauri stars (TTS in the following) is controlled by the 
interplay of different physical processes, which include vis- 
cous accretion onto the central star, photoevaporation by 
the stellar radiation field, and planet formation. If the low- 
mass star is sufficiently close to a more massive and hotter 
object, photoevaporation by this external source of high- 
energy radiation is also important because it can dominate 
the disk evolution. 

The process of photoevaporation has been extensively 
discussed (see, e.g. Hollenbach et al. PPIV; DuUemond et 
al. PPV). High-energy photons heat gas disks to temper- 
atures such that the thermal pressure exceeds the gravity 
from the central star. The disk evaporates from outside- 
in, with a mass-loss rate that decreases with time as the 
disk shrinks. Most of the well-studied Orion proplyds are 
caught in the first, short evolutionary phase, when the 
mass-loss rate is very high {Mioss > 10^^ Mo/yr, Henney 
et al, 1999). 

We report in this Letter on the first detection of a pro- 
plyd in a much later evolutionary phase, when the mass-loss 
rate is < 10~^ Mg/yr. The object (a low-mass TTS) is lo- 
cated in the star- forming region a Ori (age 2 — 3 Myr), 
which contains the massive quintuplet system a Ori. The 
brightest star of this system (spectral type 09.5) forms a 
large, low-density HII region (Habart et al., 2005) and the 
bright PDR known as the Horsehead (Abergel et al., 2003). 
The low background from the HII region and the relatively 
short distance 400 pc; Mayne & Naylor, 2008) allow 
us to detect the slow photoevaporated TTS wind and to 
measure its optical line emission spectrum. 



Send offprint requests to: erigliaco@arcetri.astro.it 

* Based on observations collected at the European Southern 
Observatory, Chile. Program 074.D-0136(A). 



2. The T Tauri star S0587 

2.1. Stellar properties 

S0587 (also identified as R053833-0236 Wolk J.S., 1996 or 
Mayrit 165257 Caballero J.A., 2008) (a2ooo=+05:38:34.04, 
'52ooo=-02:36:37.33) has been classified as M3-M4 by 
Zapatero-Osorio et al. (2002) on the basis of low-resolution 
optical spectra (Te// ^ 3300K). The optical extinction is 
negligible (Oliveira et al., 2004), and we estimate a lumi- 
nosity of 0.3±0.1 Lq, based on this spectral type and the 
V, R, I magnitudes (Wolk, 1996). The corresponding mass 
is about 0.2 M©, both from D'Antona & MazziteUi (1997) 
and Baraffe et al. (1998) evolutionary tracks. Based on its 
location on the HR diagram, S0587 has an age of '--^ 1 Myr, 
apparently slightly younger than the bulk of the a Ori stars. 
S0587 is an X-ray source, with luminosity Lx ^ 10^^ erg/s 
(Franciosini et al., 2006) 



2.2. Disk properties 

S0587 is detected by Spitzer at all IRAC wavelengths and 
at 24 fim with MIPS. It shows a relatively weak excess emis- 
sion, and is classified as an evolved disk (EV, Hernandez et 
al. 2007). We confirm this classification. The IR emission is 
reproduced well by a geometrically flat, optically thick disk 
heated by the central star, seen at an inclination of about 
40 deg. The outer disk radius is not constrained by the ex- 
isting infrared photometry, which is limited to 24 fim, as 
long as Rout ^1 — 2 AU. The model is not unique but explo- 
ration of a large number of disk models rules out optically 
thin disks, which have different Spitzer colors and lower lu- 
minosity than observed in S0587 {Ldisk/^star^ 0.07) (see 
Cieza et al., 2007). Also, there is no evidence in the SED 
of the large inner holes (few AUs) seen in transitional disks 
(Chiang & Murray-Clay, 2007 and references therein), as 
all the models that fit the data have inner radii < 0.1 AU. 
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Fig. 1. The top panel shows the spectrum in the Ha region, 
with Ha and the two [Nn]A6548.05 and [Nn]A6583.45 hnes. 
Middle panel: [Nn]A6548.05 and [Nn]A6583.45 hne profiles 
as function of the velocity shift with respect to the stel- 
lar velocity (vertical dashed line). Bottom panel: same for 
[Sn]A6716.64 and [Sn]A6730.82, as label. 
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Fig. 2. Distribution of the ra- 

tios [OI]6300/([SH]6716-I-[SII]6731) and 

[Nn]6583/([Sn]6716+fSH]6731) among the Hartigan 
et al. (1995) Taurus TTS with detected lines. The magenta 
dashed region shows the location of S0587 consider- 
ing the errors on the line equivalent widths; [01] from 
Zapatero-Osorio et al., 2002. 



2.3. Accretion rate 

S0587 shows no evidence of significant accretion onto the 
central star. An estimate of Mace can be obtained from the 
U band magnitude 18.5 ± 0.8 measured by Wolk (1996). 
For the spectral type and luminosity of S0587, this cor- 
responds to an excess U band emission of 1 mag at most. 
Using the correlation between the U band excess emission 
and the accretion luminosity established by GuUbring et 



al. (1998) for TTS, we derive Lacc^ 3 x 10"^ L© with an 
uncertainty of a factor 2, i.e., about 1% of the stellar lumi- 
nosity. Since a comparable U band excess emission may be 
due to chromospheric activity (White & Ghez, 2001), this is 
likely an upper limit to hacc- The corresponding accretion 
rate is <3 x W-^° Mg/yr. 



3. The optical spectrum 

S0587 was observed by Sacco et al. (2008) using the multi- 
object instrument FLAMES on VLT/UT2 and Ghaffe spec- 
trograph with the HR15N grating (6470-6790 A, spectral 
resolution R=17,000). The observations were obtained in 
6 separate runs of 1 hour each in October and December 
2004. Details on the observations and data reduction, in 
particular on sky subtraction, can be foimd in Sacco et 
al. (2008). The spectra of S0587 show Ha in emission, as 
well as strong forbidden lines from [NH] at 6548 and 6583 
A and [SH] at 6716 and 6731 A (Fig.[l]). Table [1] gives the 
pseudo-equivalent width (pEW) of each line (negative val- 
ues for emission lines), the FWHM (not corrected for the 
instrumental resolution of about 20 km/s), and the line lu- 
minosities. The contribution to the lines (both Ha and the 
forbidden lines) of the extended HH region due to the 09.5 
star a Ori is negligible. 

The line luminosities were computed using the observed 
R magnitude to calibrate the continuum flux and the wave- 
length dependence within the band appropriate for a 3300K 
star. The luminosity is of the order of 10^^ Lq for all the 
forbidden lines. 

Zapatero-Osorio et al. (2002) also detected Ha and for- 
bidden emission lines in two of their low-resolution spectra 
with highest quality, but not in the other, lower-quality 
spectra; when detected, the pEW of the lines are similar 
to our values. We checked for evidence of line variability in 
the Sacco et al. spectra, and found none. The line profiles 
and intensities are very similar in each run, with the pEW 
varying ^ 10%. 

In the two spectra where emission lines were detected, 
Zapatero-Osorio et al. (2002) also measured the pEW of 
the two oxygen lines at 6300 and 6364 A, which are outside 
our spectral range. The corresponding luminosity is 1.2 ± 
0.06 X 10-5 L© for the 6300 A line and 0.4 ± 0.03 x 10"^ 
Lq for the 6364 A one, similar to the [SII] and [Nil] line 
luminosity. 



4. Not an accretion-driven wind 

The forbidden lines detected in S0587 are seen in many 
TTS, and are usually interpreted in the framework of 
accretion-driven winds (e.g. Ferreira et al., 2006 and refer- 
ence therein). These centrifugally-driven MHD winds orig- 
inate in the inner disk and move fast. An estimate of the 
average mass-loss rate can be obtained from the luminos- 
ity of the forbidden lines, following Hartigan et al. (1995). 
Taking, e.g., the [SII]A6731 line, this is 



M^^nd = 3.38 X 10-*^ ( 1 



He 
He 



^6731 



150 km s 



2 X IQis cm 



lO-^L©^ 

(1) 

Me/yr. 
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The electron density rig, derived from the ratio of the two 
[SII] hnes, is ~ 1.5 x 10'^ cm~^, well below the critical den- 
sity ric, which is ~ 1.3 x 10^ cm"^. Assuming that all S 
is SII, a S abundance of 1.6 x 10~^, a typical value of 
V± ~ 150 km for the outflow velocity and that the out- 
flow fills our beam of 1 arcsec {l± ^ 6 x 10^^ cm), we derive 
a mass-loss rate of about 10~^ Mq/yt. 

This value of M^ind is high for a 0.2 M© TTS (Hartigan 
et al., 1995). It is particularly high when compared to 

Mace, as Mfjuind/ Macc^^O- In TTS, this ratio is < 1 and 
< 0.1 in the majority of cases (Hartigan et al., 1995; 
White & Hillenbrand, 2004), in agreement with the ex- 
pectations of accretion-driven jets and winds (e.g., Shu et 
al., 1994; Pudritz et al., 2007). Even taking the uncer- 
tainties on the assumptions entering into Eq.(l) into ac- 
count (see Cabrit, 2002), the very high value of the ratio 
Mujind/ Mace makes it unlikely that the forbidden lines in 
S0587 could be emitted by an accretion-driven wind. 

5. Line profiles 

Additional evidence that the S0587 forbidden lines do not 
originate in an accretion-driven wind is provided by a com- 
parison of the forbidden line luminosities and profiles with 
those of classical TTS. We first notice the strength of the 
[Nil] lines, which are weak or absent in TTS, but compa- 
rable in luminosity to the [SII] lines in S0587. Strong [Nil] 
lines are expected in highly ionized gas, as in HII regions or 
in the bow shocks at the head of jets, and are weaker than 
the other forbidden lines in the partially neutral wind/jets 
of TTS. S0587 is also unusual in having strong [SII] lines 
with respect to the [01] hnes (Fig. [J). 

Moreover, the profiles of the S0587 [SII] and [Nil] lines 
differ from the typical TTS hne profiles. The [SII] and [Nil] 
are resolved in our spectra, and they show a small shift 
of the peak to the blue (<10 — 14km s~^), a steep rise 
in the blue, and a slow decline in the red wing. Spatially 
unresolved profiles of the [SII] 6731 A line in a sample 
of accreting TTS typically show larger peak blucshifts 
(~ 80— 100km s^^), and significantly more emission in the 
blue than in the red wing (Edwards et al. 1987; Hartigan 
et al. 1995). In fact, these asymmetries are interpreted as 
evidence that the disk obscures the receding portion of the 
line emitting region. 

The high luminosity in [Nil], the low ratio [01] /[SII] 
~ 1, and the details of the line profiles of S0587 are strongly 
reminiscent of the properties of the Orion proplyds (Bally et 
al., 1998; Henney & O'Dell, 1999). The Orion proplyds are 
interpreted as the result of photoevaporation of the disk by 
the hot Trapezium stars (see, e.g. Hollenbach et al. PPIV). 
The optical forbidden lines of ionized species, such as [SII] 
and [Nil], are formed in the outer part of the outflow, which 
is ionized, heated, and shocked by the radiation and wind 
of e^C Ori (Garcia- Arredondo et al., 2001). 

6. A photoevaporating disk 

We propose that S0587 is a photoevaporating disk, which 
we detect in optical forbidden lines because it is ionized 
by the nearby (projected distance of ^--^0.35 pc) 09.5 star 
(7 Ori. Let us assume that the S0587 disk loses mass at a 
rate Mioss- Whatever the source of the photoevaporating 
photons (the central star or a Ori - see below), this wind 



Table 1. Characteristics of the emission lines. 
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originates in the outer disk, and is slow (few km s^^) and 
neutral. However, the high-energy (hv >13.6eV) photons 
from (7 Ori ionize the outflowing gas to a distance rip from 
S0587 roughly 

Here a is the recombination coefficient {a = 2.6 x 
10~^^ cm^ s~^), the mean molecular weight $i the high 
energy flux from a Ori ($i ~ 10**^ photons/s; Peimbert & 
Rayo, 1975), A is the distance between S0587 and a Ori, 
and vjp ^ 10 km s~^ the velocity at the ionization front 
(Stoerzer & Hollenbach, 1999), so that 

/ li'f, \ 2/3 / A s 2/3 

riF - 80AU f : ) f -^^) . (3) 

VlO-SMo/yr/ Vo.35pc/ ^ ' 

The ionization front moves away from S0587 and is 
closer to a Ori, if either Mioss and/or the distance between 
the two stars increases. The optical forbidden lines shown 
in Fig. [T] can only come from the ionized regions of the 
outflow; since they are detected in spectra taken with a 
beam of radius ^ 200 AU centered on S0587, it must be 
riF <^ 200 AU. This condition is verifled for Mioss ^ 10"^ 
MQ/yr; if Mioss is signiflcantly higher, rjF moves too far 
from S0587, and the outflow region within the beam is all 
neutral. 

In this model, the [SII] and [Nil] forbidden lines 
originate mostly in the ionized gas at rjF- For M;oss = 10~^ 
MQ/yr, the electron density at rjF from the continuity 

equation {ue — Mioss/{4:'Trjpfj,HViF) is ~ 2 x 10^ cm~^, 
in good agreement with the electron density derived from 
the observed ratio of the two [SII] lines 1.5 x 10^ cm~^; 
§4). The [SII] 6731 A line luminosity within the beam, 
computed assuming that all S in SII and a temperature 
^ lO'^K, is 10~^ Lq, in good agreement with the 
observations. 

What is the origin of this outflow? Both the S0587 star 
and a Ori can, in principle, cause disk photoevaporation, 
and the properties of the outflow are, qualitatively, similar. 
Let us consider first the case of disk photoevaporation due 
to the FUV (6eV< hv <13.6eV) radiation of <t Ori. Adams 
et al. (2004) have computed detailed models of the photo- 
evaporation by an external source of disks around stars of 
different masses. At a distance of 0.35 pc, the FUV radia- 
tion field from a Ori is Go ~ 10^ (in units of the Habing 
radiation field; Abergel et al., 2003). Let us assume that the 
S0587 disk has mass ^ 0.05 Mgtar when the outer radius 
is ~ 60 AU and surface density oc R~^. If we take Adams 
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et al. (2004) results for a star of central mass 0.25 Mq in 
a radiation field Go = 3000 (the exact value of Go is not 
crucial as long as it is higher than this value), the disk has a 
first, evolutionary stage characterized by a relatively large 
size and high mass-loss rates ~ few x 10~^ M0/yr. This 
phase is rapid, and the disk loses about 80% of its mass 
in ~ 3 X 10^ yr. By that time, the disk radius is Rd ^ 10 
AU, and the mass loss rate is ~ 10~^ Mq/yt. In this later 
evolutionary stage, the small, slowly evaporating disk has 
a lifetime of ~ 2 x 10^ yr, comparable to the age of S0587 
itself. 

This scenario can explain OTir data, as long as S0587 
stays in the vicinity of a Ori for ^lO*" yr. With a velocity 
dispersion of ci ^ 1 km/s in the radial direction (Sacco et 
al., 2008), the distance between a Ori and S0587 can in- 
crease to ~ Ipc in 10^ yr, i.e., about 3 times the present 
projected distance. This is still consistent with the con- 
straints we have set at the beginning of this section, al- 
though marginally so. 

However, it is possible that the photoevaporation of the 
S0587 disk is caused by the high-energy radiation from the 
central star. Recently, Gorti & HoUcnbach, (2008) model 
the case of a low mass TTS (0.3 M©), with properties 
roughly comparable to those of S0587. The trend of Mioss 
with time is similar, with an early phase when the disk is 
large and Miogs>lO~^ Mq/jy, followed by a decrease in 

the disk radius and of Mioss ■ The phase with M; 10-y 
Mo/yr, Rd ^ 30AU has a lifetime of about 5 x 10^ yr, 
which can account for the S0587 observed properties. 

In both first period of rapid evaporation, which 

produces a significant reduction of the disk size, is followed 
by a phase of slower evolution, with correspondingly lower 
values of Mioss- We think that S0587 is in this late evo- 
lutionary stage, because its disk has shrunk to 10 30 AU 
in radius and is now in the last, long phase of slow pho- 
toevaporation. S0587 shines in the optical forbidden lines 
of species, such as [SII] and [Nil], only because it is at the 
right distance from a Ori at the time of our observations. 

The exact values of Mioss, the location of the ionization 
front, the dynamics of the gas, and in particular the line 
profiles arc all the result of the interaction of several com- 
plex physical processes, that need to be explored. We note, 
for example, that the observed profiles of the optical forbid- 
den lines require a line-emitting region to be much larger 
than the disk, which does not seem to obscure the receding 
part of the outflow. Observationally, there are several pos- 
sible tests. The disk size can be constrained directly by spa- 
tially resolved observations at millimeter wavelengths, but 
also indirectly by photometry at far-infrared wavelengths. 
Also, a basic feature of externally photoevaporated disks 
is the displacement between the ionization front and the 
star; in S0587 this is predicted to be ~ 0.25 arcsec, de- 
tectable with spectroastrometric techniques or by direct 
imaging from HST. 

7. Discussion and conclusions 

The possibility that we are detecting the final, but long- 
lasting stages of the dissipation of a disk by photoevapora- 
tion is intriguing. 

As mentioned in §1, a Ori is particularly suited to this 
kind of study. First of all, the low background emission is 
very favorable for detecting objects, such as S0587, with 



a mass-loss rate at least two orders of magnitude lower 
than that of the Orion proplyds, and, more generally, for 
studying the effects of an external, high-energy radiation 
field on jets (Andrews et al., 2004) and disks. 

A second important point is the age of a Ori (2-3 Myr), 
significantly more than of Orion, such that disks have had 
enough time to evolve to the stage of low-accretion rate, low 
photoevaporation mass loss. We note that Mioss^ Macc'^ 
10"^ Me/yrin S0587, as expected at time >10'^ yr in disk 
evolution models where the photoevaporation by an exter- 
nal source is coupled to viscous evolution (Clarke, 2007). 

Some disks, such as 80587, show evidence of grain 
growth and settling in their SEDs, and it seems possible not 
only to detect old photoevaporated disks but also to relate 
their properties to those of the grains. Larger grains have 
two competing effects on the photoevaporation process: on 
the one hand, they have lower extinction in the FUV, al- 
lowing the FUV radiation to penetrate more deeply into 
the disk; on the other, the heating efficiency will be lower. 
The balance between these two competing effects needs to 
be explored in the model calculations, but can also be ad- 
dressed observationally. 

Within 0.4 pc around a Ori Hernandez et al. (2007) 
detect about 30 objects with all kinds of disks (Class 
II, Class III, and EV). The spread of SEDs probably 
reflects a varying degree of dust evolution (DuUemond & 
Dominik, 2004; D'Alessio et al., 2006). We are currently 
carrying on a comprehensive study of the accretion and 
mass-loss properties of a Ori objects to identify other 
proplyds and to address some of the issues mentioned in 
this section. 

Acknowledgments: We would like to thank S. Wolk and 
D. HoUenbach for useful discussions. We also thank the 
anonymous referee for comments that clarified the paper. 
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